Quantum critical point with competing propagating and diffusive spin excitations 
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Feedback effects due to spin fluctuation induced precursors in the fermionic quasiparticle spec- 
trum are taken into account in the description of a quantum critical point of itinerant spin systems. 
A correlation length dependent spin damping occurs, leading to a dynamical scaling with z « 1 
which non-trivially competes with the conventional spin wave behavior. We obtain, within a one 
loop renormalization group approach, a quantitative explanation for the scaling behavior seen in 
underdoped cuprate superconductors. 



Nuclear magnetic resonance (NMR) and inelastic neu- 
tron scattering (INS) experiments in magnetically under- 
doped cuprates [l]-[|] indicate that the spin excitations 
of these systems are close to a zero temperature quan- 
tum critical point (QCP) These measurements 
give strong indications for a dynamical scaling exponent 
z « 1, typical for systems characterized by propagating 
antiferromagnetic (AF) spin modes, but different from 
that expected for an itinerant AF system with z = 2 
H^l]] ■ On the other hand, the low frequency behavior 
of the imaginary part of the dynamical spin susceptibil- 
ity is characterized by over-damped spin-excitations, i.e. 
Imx(q, u>) oc ui + 0(uj 3 ), as expected for itinerant AF. 
Thus, while the spin dynamics is clearly over-damped, 
the scaling behavior corresponds to that of a system with- 
out damping, causing a conceptual problem for the de- 
scription of the related QCP. 

In the present paper we investigate the behavior in the 
vicinity of a QCP characterized by competing propagat- 
ing and diffusive spin excitations. We explicitly take into 
account that for a system characterized by strong interac- 
tions between collective spin modes and fermionic quasi- 
particles the latter are strongly affected by the critical 
mode. This, in turn, leads to a feedback in the spin dy- 
namics [ ^0|JTl| ] , causing zwl dynamical scaling behavior 
for an itinerant AF. Our theory offers a quantitative ex- 
planation for the temperature and frequency dependence 
of INS and NMR data and confirms the phenomenologi- 
cal description of the spin dynamics of Refs. . 

We start from an interacting Fermi system charac- 
terized by some unperturbed band-part and an interac- 
tion term, X) q /q s q ' s -q; with fermionic spins, s q = 

\ Skcrcr' c k+qo-°'ffff' c kff'- By introducing a collective spin- 
1 Bose field, S(q), via the Hubbard-Stratonovich trans- 
formation, one can integrate out the fcrmions and expand 
up to forth order in S(q) fl2|| ; the resulting effective ac- 
tion of the collective spin degrees of freedom is : 
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S=\ fdqx^iq) S(q)-SHz) 



I <%•••/ dqi S qi+q2+q3+q4 
xS(q 1 )-S(q 2 )S(q 3 )-S(q 4 ), (1) 

characterized by a bare spin propagator, Xo{q)i an d cou- 
pling constant, u > 0. In Eq. [I] the (d + I)-dimensional 
vector q = (q, iu) n ) consists of the d-dimensional mo- 
mentum vector q and the bosonic Matsubara frequency 
io n = 2nnT at temperature, T. We write J dq ■ ■ ■ = 

T J2 n Jj qJ <A iCf ' ■ • and <W = r-'W^Hq + q'), 

where A is the upper momentum cut off. A diagrammatic 
representation of this action is given in Fig. [j](a). 

Usually, Xo(q) an d u are calculated within a weak cou- 
pling approach in which the spin damping due to particle 
hole excitations is identical to that of a noninteracting 
electron gas (see Fig. 0(b)). It is however to be 

expected that the fermionic quasiparticles are strongly 
affected by their own spin excitations once the system 
gets close to a magnetic instability JTT| , ^3| -^5| . It is 
therefore interesting to take into account these modi- 
fications of the quasiparticle spectrum in an investiga- 
tion of the critical behavior. Of particular importance 
are quasiparticles close to the magnetic BZ-boundary, 
which, on the one hand, determine the spin damping, 
7 oc Imx _1 (Q, w)/w| w ^0j and, on the other hand, are 
mostly affected by the proximity to the ordered state. 

Magnetic precursors in the quasiparticle spectrum have 
recently been investigated in several theories jll],[l3] |l5| . 
Within a quasistatic theory ftijj , applicable at higher T, 
7 was calculated to all orders in the perturbation the- 
ory. Due to precursors of a spin density wave (SDW) 
gap the low energy spectral weight is reduced and the ver- 
tex function enhanced and 7 was found to change from 
being a constant to a ^-dependent function: 7 cx 
with ip fa 1, once £ is larger than the electronic length 
scale, £ = 2w//q. In cuprates £ ~ 2 due to the large 
coupling constant /q. For low temperatures the sit- 
uation is more complicated p6| , p7t . Nevertheless, the 
leading contribution to the spin damping can be esti- 
mated as 7 oc /qT /v 2 , with spin fermion vertex func- 
tion at the Fermi surface (FS) and for momentum trans- 
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for Q w (tt, 7r), T sf , and Fermi velocity, v. At the critical 
point the vertex function for scattering processes includ- 
ing a Goldstone mode is expected to vanish on a surface 
defined by momentum and energy conservation of the in- 
volved quasiparticles pl| . In the case of large correlation 
length and large interaction /q, a FS evolution towards 
an SDW like behavior occurs and the above principle 
causes the vertex function to vanish like T s f oc [ fL4| , 
leading to 7 oc Thus, we expect in various regimes 
a ^-dependent reduction of the spin damping due to the 
proximity to an ordered state. 
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FIG. 1. (a) Diagrammatic representation of the effective 
action with spin-propagator and spin spin interaction term. 
The dashed lines correspond to collective spin degrees, the 
solid lines to the renormalized fermionic Green's function, (b) 
Typical weak coupling approximation for Xo{l)- Dotted lines 
refer to the fermion-fermion interaction, / q . (c) Some of the 
diagrams causing feedback effects of the quasiparticle excita- 
tions. 

In the present communication we take these feedback 
effects into account and determine £(T) using a renormal- 
ization group approach. The propagator of the collective 
spin modes is assumed to be: 



xo(q,«w„) 



£0 +q 2 + £ ¥ >«l/c+K/c) 5 



(2) 



with bare correlation length £01 spin damping, 7 = 
c _1 ^o~ ¥ ' an( i spin wave velocity c. Momenta are measured 
relative to the ordering vector. For c — > 00, i.e. without 
spin damping, the problem is similar to that investigated 
by Chakravarty et al. Q and Chubukov et al. ||. The 
alternative limit ip = and c — > 00 was discussed by 
Hertz Jj| and Millis (7]]. The case ip = but c finite was 
discussed by Sachdev et al. J|], who found a z = 1 to 
z = 2 crossover for decreasing temperature, in contrast 
to the experimental finding in the cuprates of a z = 2 to 
z = 1 crossover for decreasing T @,|,|. 

We use a Kadanoff- Wilson momentum shell renor- 
malization group approach; integrating out states with 
momenta between Ae~ l and A, including the rescaling 
T(Z) = e zl T to reach self-similarity. Up to one loop the 
resulting flow equations are: 

dT(l) 



dl 

du(l) 

dl 
dc(l) 

dl 
dc(l) 

dl 



= 2e o - 2 (0 + 4(A + 2H>(£,T(Z)) (4) 

= €u(l)-4(N + 8)u 2 ^(tT(l)) (5) 

= (z + <p- 2)c(0 (6) 

= {z- l)c(i) , (7) 



dl 



zT{l) 



with e = 4 — (d + z). N is the number of vector compo- 
nents of S(q). The functions <f> and "J can be determined 
along the same lines as in Ref. Q . Their explicit depen- 
dence on the renormalized correlation length, £, results 
from a replacement of £0 by £ in diagrams like those in 
Fig. |l](c) which determine the spin damping, necessary 
to reach self consistency |l9|]. From Eqn. ^ and we find 
that for ip = 1 (i.e. z = 1) the velocities c and c do 
not renormalize in the one loop approximation. In the 
case ip < 1 one starts from a primary scaling behavior 
with z = 1 which eventually crosses over to z = 2 — ip 
||. In what follows we assume ip = 1 since quantita- 
tively the latter crossover will barely change our results; 
an extension to arbitrary ip is straightforward. More- 
over, small deviations from z = 1 seem to be beyond 
the accuracy of the current experiments. For e > 
the flow equations, Eqn. || - 0, are characterized by a 
zero temperature fixed point (£o) -2 = ~W^k e ^ 2 anc ^ 
u* = T(d/2)2 d - 1 -K d ^e/((N + 8)cA d ~ 3 , which is unsta- 
ble since T is a relevant field. By integrating the flow 
equations we find that the critical behavior is sufficiently 
characterized by the single effective coupling constant 
g = 2(1 - d)(N + 2)u^. This is also expected since the 
critical behavior should be similar to that of the quan- 
tum nonlinear sigma model which is fully characterized 
by g and T (|||. For d — 2, the fixed point value of g 
is given by g* = 47r/(cA), which is unchanged from the 
result without damping JsjJ^] (the limit c — > 00). Never- 
theless, the critical behavior in the vicinity of this QCP 
depends strongly on the ratio, T = c/c. In the limit 
c — > 00 the QCP moves towards g* — 87r/A w with upper 
frequency cut off A u . The flow equations, Eqn. || - [7], 
can be integrated and the correlation length follows from 
£ = e £(<?(/), T(£)), where we determine the correlation 
length at the matching point ||, I = I, with T(l) = cA/4 
and g(T) = 2g* from a l/N expansion as used in Ref. ||. 
Note, since g(l) depends on £, this procedure leads to a 
self consistency condition for £. 

The QCP, <?*, separates in the usual sense [|| a renor- 
malized classical (RC) regime with exponentially growing 
correlation length, £(T) = c/(2T) exp (2ttcA/T(^- - 
for g(l = 0) < 5*, from a quantum disordered (QD) state 
with a finite zero temperature correlation length, £(0), 
(see Fig. |). We find for T > 1 that £(0)/log(£(0)A) = 
J^A (5 — 9*)~ 1 i w bich is enhanced compared to the situ- 
(3) ation without damping: £(0) = fx (.9 — with z/ = 1 
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up to one loop. The results of a numerical evaluation of 
^(0) _1 are shown in the inset of Fig. ||. 



the momentum averaged susceptibility J d 2 qlmx q (w) is 
found and the low frequency slope of X q =o(w) behaves 
like ImXq=o('^)/w|( 1 ;— >o = a£, 5 /c with S — 3. 
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FIG. 2. Phase diagram of an itinerant AF with 2 = 1 scal- 
ing behavior for T = 4. RC is the renormalized classical 
regime with exponentially large f. QC is the quantum crit- 
ical regime and QD the quantum disordered regime. In the 
mean held regime, T > T cr , with £ < £ w 2 no feedback 
effect due to changed quasiparticle behavior occurs, leading 
to z = 2. The inset shows the inverse zero temperature cor- 
relation length in the QD regime for different V. 

Modified scaling behavior can in similar fashion be 
found as function of temperature. In Fig. [| we show 
our results for the T-dependent correlation length for 
g/g* = 1.18 and different values of T. We find up to 
logarithmic corrections, £(T) -2 « £(0) -2 + bT 2 (see in- 
set of Fig. ||(b)). In contrast, in the case of propagat- 
ing spin excitations a sharp transition to an expnentially 
weak temperature dependent £ occurs for T rj A with 
spin wave gap A = c£(0) -1 . Thus, due to spin damping 
the spin wave gap is filled with low energy states and the 
correlation length continues to grow even for very low 
temperatures. The crossover between the quantum crit- 
ical (QC) and QD regime is very gradual. In Fig. || the 
crossover line was determined by the temperature where 
d 2 ^ 1 /dT 2 becomes small compared to its low T value, 
i.e. where £ -1 starts to grow linearly with T. Finally, 
considering higher temperatures, we have to take into 
account that 7 oc occurs only for £ > £ Q « 2 [ [Ti"|| . 
For £ < £ , 7 = const, leading to a z = 2 behavior at a 
characteristic temperature T cr (see Fig. ||). 

Another characteristic phenomenon caused by the 
proximity to a QCP is the scaling behavior of the fre- 
quency dependence of Imx(q, ui). For any point in the 
(T, g) phase diagram (except g < g* and T = 0) Imx q (a;) 
increases linearly with oj due to excitations in the par- 
ticle hole continuum. In the RC regime however, the 
spin damping is exponentially suppressed and the spin 
dynamics is indistinguishable from a system with purely 
propagating spin waves. On the other hand, in the QC 
regime ui/T scaling behavior of £ ImXq=o(w) and of 
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FIG. 3. Temperature dependence of the correlation 
length for g — 1.18g* and different T in comparison 
with experimental data for Lai.86Sro.i4Cu04 from Ref.[3] at 
io — 2.5 (3.5) meV indicated by diamonds (triangles) and the 
phenomenological theory of Ref. [2]. The inset shows £(T) -2 
versus T 2 . 

In applying these results to cuprate superconductors, 
we assume that doped systems without long range or- 
der are located in the QC and QD regime. Since g oc 
l/(S(r) 2 }, as follows from a 1/N expansion, doping re- 
duces the effective moment, causing the coupling con- 
stant to grow until it exceeds g*. For the cuprate ma- 
terial Lai.86Sro.i4Cu04 the parameters c « 50meV and 
c « 220meV are reasonably well known from NMR and 
INS experiments (l]-g] and it follows r w 4. Thus, it 
suffices for a quantitative understanding of the INS data 
of Ref H to determine the ratio g/g* from the correla- 
tion length at a given temperature. Once this number 
is determined one has a complete description of the T- 
and w-dependence of the spin fluctuation spectrum. As 
shown in Fig. [| a value g = 1.18g* gives a reasonable 
agreement not only for the value of £ at low T but also 
for the whole temperature regime up to 300 K. Using 
this value for the coupling constant we can determine the 
frequency and temperature dependence of the dynami- 
cal spin susceptibility. The results for Imx q= o(aj) with 
a = 26 eV -1 , shown in Fig. ^, are in remarkable agree- 
ment with the results of Aeppli et al. m , who also find 
Imx q =o( w )/ w U^o x £ S with S = 3±0.3. This agreement 
between theory and experiment is a direct consequence of 
the fact that the experimental data show a uj/T scaling 
behavior of ^ 2 Imx q= o(w) which is a strong indication 
that the system under consideration is indeed close to 
a QCP. In order to demonstrate that our results agree 
quantitatively with the phenomenological description of 
the 2 = 1 pseudo-scaling regime of Ref. [Q , we also show 
in Fig. U the correlation length, £(T), as obtained from 
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FIG. 4. Frequency dependence of the dynamical suscep- 
tibility at the peak maximum for T = 4 and g = 1.18<?* in 
comparison with INS experiments from Ref. [3]. 

In conclusion, based on previous calcula- 
tions [ pd]p^ -[l5|], we have argued that modifications of 
the particle hole excitation spectrum can change the dy- 
namical scaling behavior of itinerant AF systems close 
to a QCP by affecting the spin damping. The resulting 
z « 1 scaling causes a T-dependence of the AF correla- 
tion length which is completely different from the usual 
z = 2 case. This new T-dependence of the dynami- 
cal spin susceptibility is in remarkable agreement with 
NMR and INS experiments and gives an explanation for 
the crossover scenario of Refs. 0,^]. The z = 1 scaling 
and the position of the QCP are the same as those for 
propagating spin modes in insulating AF. However, due 
to damping, the correlation length is enhanced and the 
T-dependence of £ is changed. 

We expect similar behavior in the quantum phase tran- 
sitions of other itinerant systems if the strong interaction 
between quasiparticles and collective modes changes the 
dynamics of the collective mode under consideration. Ex- 
amples are one and two dimensional charge density wave 
systems, superconductors or saturated ferromagnets. In 
higher dimensions or for systems with weak quasiparticle- 
collective mode coupling precursor phenomena are only 
expected once £ diverges |lq ] and will play no role in the 
quantum disordered regime. 

Finally, we note that the theory presented in this pa- 
per is not a theory for the crossover to a strong pseudogap 
state found in many cuprate superconductors at low T, 
where a decoupling of the T-dependence of £ and 7 has 
been found UH. This might be due to the interference 
with excitations different from the spin fluctuation chan- 
nel discussed here [|!7|j2lf| . Nevertheless, we expect the 
proximity to a QCP to be essential for the appearance of 
the strong pseudogap behavior. 
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